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Abstract. We discuss the stick-slip motion of an elastic block sliding along a rigid substrate. We argue 
that for a given external shear stress this system shows a discontinuous nonequilibrium transition from a 
uniform stick state to uniform sliding at some critical stress which is nothing but the Griffith threshold 
for crack propagation. An inhomogeneous mode of sliding occurs, when the driving velocity is prescribed 
instead of the external stress. A transition to homogeneous sliding occurs at a critical velocity, which is 
related to the critical stress. We solve the elastic problem for a steady-state motion of a periodic stick-slip 
pattern and derive equations of motion for the tip and resticking end of the slip pulses. In the slip regions 
we use the linear viscous friction law and do not assume any intrinsic instabilities even at small sliding 
velocities. We find that, as in many other pattern forming system, the steady-state analysis itself does not 
select uniquely all the internal parameters of the pattern, especially the primary wavelength. Using some 
plausible analogy to first order phase transitions we discuss a "soft" selection mechanism. This allows to 
estimate internal parameters such as crack velocities, primary wavelength and relative fraction of the slip 
phase as function of the driving velocity. The relevance of our results to recent experiments is discussed. 

PACS. 62.20.Mk Fatigue, brittleness, fracture, and cracks - 46.50.+a Fracture mechanics, fatigue and 
cracks - 46.55.+d Tribology and mechanical contacts - 62.20.Qp Tribology and hardness 



1 Introduction 

The nature of sliding friction and especially of inhomoge- 
neous modes of sliding is a fundamental physical problem 
of prime practical importance 1 . Studying the friction 
of rubber along a smooth glass substrate, Schallamach |2] 
observed an inhomogeneous mode of sliding. Under some 
conditions homogeneous sliding becomes unstable and a 
quasiperiodic pattern of detached zones is formed. Re- 
cent experimental observations of Rubio and Galeano [3], 
Baumberger, Caroli, and Ronsin on the frictional 

motion of sheared gels sliding along a glass surface also 
indicate the existence of self-healing pulses and inhomoge- 
neous modes of sliding. A regime of periodic stick-slip has 
been observed in a limited range of small shear rates 01 
[5] . It bifurcates towards stationary sliding at some critical 
driving velocity. The slip pulses traverse the sample with 
a velocity much larger than the driving velocity but still 
much smaller than the speed of sound. The importance 
of the slow crack-like fronts for the onset of the frictional 
slip was also stressed in [H]. 

Slip pulses in gels seem to be very different from Schal- 
lamach waves and "brittle" pulses studied by Gerde and 
Marder [7] since no observable opening occurs at the in- 

a Present address: Institut fur Theoretische Physik, Univer- 
sitat zu Koln, 50937 Koln, Germany 



terface. In this respect, they are more comparable to self- 
healing cracks, suggested by Heaton 8 in the context of 
seismic events. 

Stick-slip motion with the transition to sliding above 
some critical velocity was also observed in studies of fric- 
tion in an ultrathin layer of lubricant between two atom- 
ically flat surfaces 9 . This behavior was attributed to 
the confinement-induced freezing of the lubricant and its 
melting due to shear stress. 

Recent investigations (see, for example, ^U] and refer- 
ences therein) point towards an essential importance of the 
underlying friction law in the slip state. It has been proved 
that the simple Coulomb friction leads to the so-called "ill- 
posedness" of the linear stability problem for small inho- 
mogeneous perturbations of the stress and strain fields in 
a sliding mode (TJ3|- Moreover, Caroli JI] has shown that 
the existence of slow, periodic slip pulses is incompatible 
with the Coulomb friction law. 

In our previous paper JT2] we developed a conceptually 
simple model which is compatible with an inhomogeneous 
mode of sliding, the existence of slowly moving slip pulses 
and the "shear melting" phenomenon. We discussed the 
propagation of a shear crack (Mode II crack) along the in- 
terface between two dissimilar materials. The crack edge 
separates two states of the interface, "stick" and "slip" . 
We assumed that the interface is flat with a strong ad- 
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hesion contact. In the presence of roughness, the assump- 
tion of strong adhesion and full contact at the interface 
is presumably only reasonable for "soft" materials with a 
relatively small shear modulus. Gels are clearly materials 
of this type. 

In the slip region we assumed a simple linear viscous 
friction law, namely, that the shear stress is proportional 
to the sliding velocity. This, from the theoretical point of 
view strongly motivated law, is usually not discussed in 
the literature since it does not lead to the so-called static 
friction phenomenon observed experimentally. However, 
in our description static friction appears in a natural way 
as the usual Griffith threshold for crack propagation. The 
important point is that, before the system goes into a 
sliding mode, a shear crack should traverse the sample. 
This requires a finite shear stress since the stick state of 
the interface is energetically more favorable. 

We note that the above mentioned interface proper- 
ties in our model are fully coupled to the bulk elasticity 
of a gel leading to a rich spectrum of phenomena where 
two intriguing problems, crack propagation and interface 
friction, come together. 

In the present paper we attempt to find a solution of 
this model which represents a periodic pattern of stick-slip 
motion. The paper is organized as follows: 

In the next section we formulate the problem of the 
periodic pattern which appears in an elastic solid block 
sliding on a fiat rigid substrate (see Fig. QJ. We formu- 
late the boundary conditions for this elastic problem and 
introduce the thermodynamics of two states of the inter- 
face: stick and slip. Each slip pulse can be considered as 
a mode II crack which is bounded by two crack edges: 
the pulse tip and the resticking end. In order to formu- 
late equations of motion for such cracks we have to solve 
the elastic problem and calculate the energy flux into the 
crack edges. 

In section [3] we solve the elastic problem in two lim- 
iting cases: the height of the block is much larger than 
the lateral scale of the pattern and the opposite limit. 
This section is rather technical but eventually provides 
expressions for the energy flux into the crack tips. It turns 
out that, as in many pattern forming systems (for exam- 
ple, directional solidification or eutectic growth the 
steady-state description itself does not allow to select the 
primary wavelength of the pattern. We will see that the 
degeneracy of our system is even higher than in crystal 
growth problems. 

In section 01 we give some plausible arguments based 
on the analogy to the first order phase transitions which 
lead to the "soft" selection of the periodic pattern. We 
predict a regime of stick-slip motion in a limited range of 
small shearing rates and a transition towards stationary 
sliding at some critical driving velocity. In section [S] we 
discuss our results and their possible relevance to existing 
experimental observations. Technical details are presented 
in the Appendix. 



2 Formulation of the problem 

Consider an elastic solid sliding on a flat rigid substrate. 
Assume that the elastic solid occupies the space < y < 
H, and let (x, y, z) be a coordinate system with the plane 
y = corresponding to the surface of the solid, see Fig. ^ 
We discuss the plane strain situation with u z = 0, where 
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Fig. 1. (top) an elastic body sliding on a rigid substrate; (bot- 
tom) motion of a single crack tip. 

u is the displacement vector and static elasticity since all 
motions are expected to be slow compared to the sound 
velocities. 

Let us discuss the boundary conditions. First of all 
we assume that displacements at the top boundary of the 
block are the same as for the top of the rigid plate: u y — 
and u x — vt + const are uniform for y = H. Here v 
is the driving velocity of the upper rigid plate relative 
to the bottom rigid substrate; t is time. We assume also 
that the bottom interface can be in two states: "stick" 
and "slip". The boundaries between these two states are 
described by the crack edges which move with a velocity 
c in the x direction. In the stick regions u v — and the 
sliding velocity also vanishes, ii x = for y = 0. In the 
stationary situation, when all stick-slip boundaries drift 
with the same velocity c, displacement derivatives depend 
on x and time t only in the combination x — ct. Therefore, 
in the stick regions there is a constant deformation u xx = 
Vo = v/c. 

In the slip regions we assume that the two solids 
are always in contact, u y = for y = 0, while we al- 
low for a finite relative sliding velocity u x . This sliding 
velocity leads to frictional shear stress at the interface, 
°~xy = o~ xy {u x (t,x)}. In the following, a simple linear vis- 
cous friction law will be considered, a xy = (3u Xl with (3 
being the friction coefficient. 

Let us turn now to thermodynamical aspects of the 
problem Since the adhesion contact in the stick region is 
stronger, it is reasonable to assume that the interface en- 
ergy in the stick phase is smaller than the interface energy 
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in the slip phase. We denote this energy difference by 7. 
It is clear that without external loading the stick phase is 
energetically favorable and a finite shear stress is required 
to get the interface into the slip state. 

We first discuss the situation where the external shear 
stress <j xy is given (not the driving velocity). Then we 
have two homogeneous solutions to our problem: i) uni- 
form stick with pure elastic response and no drift velocity 
and ii) uniform sliding with drift velocity 

v = Sa xy /n, (1) 

where [i is the shear modulus and S — fj,/ (3 is the velocity 
scale given by the friction law. We note that this homo- 
geneous sliding mode is linearly stable for any velocity 
with respect to small inhomogeneous perturbations of the 
stress and strain fields. In this respect the viscous friction 
law is very different from Coulomb friction which leads to 
a linear instability and ill-posedness of the problem as has 
been intensively discussed in the literature 

In our previous paper we considered a body uniformly 
strained with u^ y at x — ► 00 and relaxed at x —> —00 due 
to the presence of the slip at x < (a singular solution of 
the static theory of elasticity similar to a crack solution, 
see Fig. ^| . The condition for the boundary between semi- 
infinite slip and stick regions to propagate in the positive 
x-direction is that the energy release because of the stress 
relaxation is larger than the energy flux to the crack tip 
(which is proportional to the surface energy): 

4 = 2/i(u~) 3 JJ/ 7 >l, (2) 

where H is the height of the elastic block. Otherwise, the 
crack would propagate with negative velocity and the stick 
phase would be restored. Equilibrium corresponds to 

<r~ = a c = V^l/H. (3) 

The condition A = 1 or a xy — a c is nothing but the 
usual Griffith threshold for crack propagation. On the 
other hand, in the context of the friction problem, this 
condition may be interpreted as a static friction thresh- 
old: a finite shear loading is required to get the system into 
the sliding mode. The energy released by the creation of 
a single slip of length I » Hj \J\ — 2v in the body, uni- 
formly stressed with u xy , is (2/z (u xy ) 2 H — Uniform 
stick in the long sample is stable if the applied shear stress 
does not exceed tr c , otherwise stick is only metastable. We 
note that for soft materials such as gels the characteris- 
tic length H/y/1 — 2v can be quite large. For length scale, 
shorter than this length, er c increases by factor of 2 com- 
pared to Eq. (PJ). We will return to this point later when 
accurately solving the elastic problem in Appendix IU1 

On the other hand, the homogeneous sliding mode may 
be unstable against a resticking pulse (nonlinear "healing 
instability") if the corresponding value of a xy < a c . Since 
in this case the shear stress is related to the steady-state 
sliding velocity by Eq. QJ, we presume that the homoge- 
neous sliding is stable against the healing instability only 
above the critical sliding velocity v c , 

v c = Sa c j[i. (4) 



Thus, for given a we discuss a dynamical first order phase 
transition between homogeneous stick and slip regimes 
with a c being the transition point (see Fig. Another 
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Fig. 2. Two homogeneous regimes, stick and slip; a c is the 
transition point and also the average stress in stick-slip regime. 

example of such a non-equilibrium transition is a dielec- 
tric breakdown above some critical value of the applied 
electric field. The chain reaction of the ionization process 
takes place forming a conducting state instead of a non- 
conducting one realized below the critical value of the ap- 
plied electrical fields. We can also mention a phenomenon 
of shear banding in complex fluids. Note that these tran- 
sitions are not real thermodynamical transitions since the 
slip phase (or conducting state) is not a real thermody- 
namical phase and dissipation due to friction takes place 
in this regime. Therefore, our arguments that a c is a tran- 
sition point are not so strict as in usual thermodynamics. 
However, we believe that this picture is plausible and also 
in the spirit of "shear melting" in ultrathin lubricant lay- 
ers, discussed in .14 . In our case, however, the surface 
properties are strongly coupled to the bulk elasticity and 
the critical stress a c depends not only on the interface 
energy difference, 7, but also on geometry. This is the 
usual situation in fracture problems and, as we have al- 
ready mentioned, a c is just the Griffith threshold for crack 
propagation in this context. 

If one starts from the stick phase and increases the 
stress slightly above the critical value, the transition to 
the slip regime goes via a nucleation process of slip re- 
gions and their growth. Eventually annihilation of the 
boundaries brings the system into the homogeneous slid- 
ing mode. If the sample is not too long, propagation of 
just one crack may be sufficient to bring the system into 
the sliding mode. The velocity c of such an isolated crack 
has been calculated in our previous paper |12| : 

C ~ wHH/a) b - [b) 

Here w is a number related to the Poisson ratio v: 

3-4^ 

and a is the small length scale cutoff (we will return to 
this point later). 
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If one starts from the sliding mode and reduces the 
stress below the critical value, the transition to the stick 
phase is via nucleation and growth of the resticking re- 
gions. 

We note that, for a given shear stress, this picture does 
not allow for a steady periodic stick-slip motion. The tran- 
sition between homogeneous stick and slip states corre- 
sponds to the critical value of the stress but, of course, 
there may be some hysteresis around the transition point, 
c c . 

The situation changes drastically when, instead of the 
shear stress, the driving velocity is prescribed. Actually, 
most of the experimental setups fix the driving velocity 
v. In this case the uniform stick state is no longer pos- 
sible. On the other hand, if the driving velocity is below 
its critical value, given by Eq. J3J|, uniform sliding may be 
unstable against resticking pulses (nonlinear "healing in- 
stability" ) . Then we expect the system to exhibit an inho- 
mogeneous sliding mode, namely periodic (quasiperiodic 
or or even chaotic) stick-slip. 

At this point the analogy to the first order phase tran- 
sitions can be made even deeper: In thermodynamics, if 
the volume but not the pressure of the system is prescribed 
for a given temperature, phase coexistence occurs in a 
range of volumes. If the volume is changed, the pressure 
remains constant according to the equilibrium phase di- 
agram and only the relative fraction of the two phases 
changes. If the volume reaches its critical value, one phase 
disappears and phase equilibrium is not possible for fur- 
ther volume changes. Of course, the driving velocity, which 
is a dynamical variable, cannot be mapped directly onto 
the thermodynamical variable, the volume, as well as the 
hydrostatic pressure is not the same as the shear stress. 
However, we use this plausible analogy and assume that 
stick-slip motion (coexistence of two phases) will occur 
in a range of driving velocities, v < v c , with the average 
shear stress a ~ cr c , as shown in Fig. [3 

Our main goal is to discuss this inhomogeneous mode 
of sliding. In this paper we analyze this problem only for 
the steady-state periodic stick-slip motion schematically 
presented in Fig. The first step is to solve the elastic 
problem for this geometry and then to calculate the en- 
ergy flux into the crack tip and the resticking end which 
allows to formulate equations of motion for these bound- 
aries. This rather technical issue is discussed in detail in 
the next section. But before some remarks are in order. 

The solution of the elastic problems defines a family 
of the displacement fields and slip velocity fields, each of 
which is labeled by five parameters: driving velocity v, av- 
erage shear stress a, crack velocity c, two length scales of 
the slip and stick regions, X s i and X s t (or alternatively, the 
wavelength of the pattern A = X s i + X s t and the relative 
fraction of the slip rj = X s i/X). The driving velocity v is 
imposed. This defines a problem of a dynamical selection, 
namely, if a sliding pattern exists, are c, X s i, X s t and a 
uniquely defined when v is fixed, or not? To find these 
four parameters we expect to have only two equations of 
motion for the crack tip and resticking end. Thus, our 
steady-state periodic problem has a high degeneracy. One 



needs two additional equations to predict a uniquely se- 
lected pattern. Using the analogy to the first order phase 
transitions we have already removed one degree of free- 
dom, selecting the average shear stress a to be close to 
its critical value a c . However, as in many pattern form- 
ing systems the selection of the primary wavelength A is 
a highly nontrivial issue. In directional solidification and 
eutectic growth (see, for example, ^3]) it is believed that 
there is only a "soft" selection mechanism which depends 
on the initial conditions, the history of the process, the 
level and origin of fluctuations, details of the nucleation 
and annihilation processes. A detailed discussion of this 
issue is far beyond the scope of this paper, but we use an 
additional assumption which leads to a "soft" selection of 
the wavelength: we assume that the smallest length scale 
among X s i and X st is of the order of the sample height H. 



3 Solution of the elastic problem and 
equations of motion for cracks. 

In this Section we solve the problem of linear elasticity 
with the boundary conditions specified above. We find the 
expressions for the energy fluxes into the crack edges in 
terms of the imposed boundary conditions (driving ve- 
locity v, average shear stress r* = er, crack tip velocity 
c) and geometrical parameters (sample height H, pattern 
wavelength A = 27r/fc and slip fraction A s ;/A = rf). We 
restrict our attention to the very slow modes assuming 
that crack velocities are much smaller than the shear wave 
speed. Stationarity of the problem shows in the fact that 
the spatial coordinate x enters only in combination with 
time: x — ct. To simplify formulas, we will usually omit the 
time dependence in the combination. 

In Appendix lAl we consider stationary solutions of the 
linear elasticity theory in the given geometry and with 
boundary conditions u y \ y= o = u y \ y= u = d x u x \ y= H = 0. 
We treat a so called mixed problem of elasticity: there are 
parts of the bottom surface (slip regions) where only a 
linear relation between shear strain and stress is given. In 
the particular case of zero friction, shear stress is zero. 

Using the solutions, we can express any relevant quan- 
tity (strain or stress components) in terms of the displace- 
ment u x \ y= o, or equivalently in terms of the dimensionless 
sliding velocity V = u x /c = Vq — d x u x \ y= Q. The relations 
are generally nonlocal, leading to an integral equation for 
the function V. 



3.1 Small wavelengths, A << H. 

Let us consider the case where the height is the largest 
geometrical parameter. Then, as in the case of infinite 
height, we obtain the following integral equation 

d 

t(0 = r* + P.V. [ dC l -^m , (6) 

71-7x7 J £ - C 

-d 
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where r is the shear stress at y = and 

3-4i/ 



2(1-1/) 



Here we made the variable transformation £ = tan(fca;/2), 
which projects a period — 7r < fcx < 7r onto the infinite 
interval — oo < £ < oo, and defined a new function = 
V^(2 arctan£)/(l + £ 2 ) which is nonzero in the slip region 
|£| < d = tan(cv/2). The average dimensionless velocity is 
given by Vq = v/c. 

First we consider the solutions of the equation in the 
absence of friction in the slip region k\x\ < a, which means 
t(£) = for |£| < d (rj = a/ir = X s i/X is the fraction of the 
slip) . The solutions were discussed in Ref. ^I] (even with a 
more complicated form of the viscoelastic bulk response): 



*o(0 



VTTd 2 wt*^ + 2fiV 



2^ (i+e)V(d 2 -e) 

The solution (JJJ can be rewritten for the strain: 

wt* tan(fca;/2) + 2^iV 

o x u x \y = o = Vo 



(7) 



2/icos(a/2)^/tan 2 (a/2) - tan 2 (fcx/2) 

(8) 

in the slip region | tan(fc.x/2)| < tan(a/2), and d x u x \ y —o = 
Vq in the stick region. The singularities in the points \x\ = 
a/k are similar to those in a crack solution in 2d elasticity 
theory. Generally, there is a macroscopic elastic energy 
flow to a point with the square-root singularity. 

The energy flow is determined by the angular depen- 
dence of the solution in the vicinity of the point. Thus, we 
can apply the result for the semi-infinite slip J2| to find 
the energy flows to the head tip J\ and from the rear tip 
Ji of the slip: 



Jl,2 = 



7rc(2/iVo ± wt* tan(a/2)) 2 
Awfik tan(a/2) 



(9) 



We see that J2 vanishes only if 2/iVb — wt* tan(a/2) = 0, 
i.e. when the sliding velocity tends to zero near the rear 
tip. This special solution was discussed in the literature 
(see, for example, Refs. [Tl"lll5lllfi| and references therein). 

The energy flowing to the tips can be reversibly con- 
verted into the interface energy and can be dissipated: 

Ji = 7 c(l + c/co), J 2 = 7c(l - c/c ), (10) 

where Cq describes tip dissipation. Let us assume that the 
tip velocity is so small that the dissipation can neglected. 
The resulting conditions, 

Ji = J2 = 7 C 



are fulfilled if either Vq — y / wjkt&n(a/2)/ /j, and r* = 0, 
or t* = 2\J^j\ikj (w tan(a/2)) and Vq — 0. The latter case 
corresponds to the periodic array of slip regions moving 
with velocity c; in the co-moving frame the strain distri- 
bution is the same as in the statics. This solution is unsta- 
ble. For example, for the fixed wave number k and stress 



smaller than its equilibrium value, a will decrease until the 
stick conditions hold everywhere. The former solution has 
the feature that the sliding preserves its direction (V > 0) 
in the whole slip region. 

The case with non-zero friction can be treated as 
above. Now, in the left-hand side of Eq. ©, in the slip 
region |£| < d, we have r(£) = /Ltc(l + £ 2 )<Z>(£) / S . The so- 
lution can be guessed from that for the semi-infinite slip 
Q2|: 



*(0 = 



(l + £ 2 )(£ + dj^+^fd - (1 " e)/2 



(11) 



where 61, 62 and e are given by the following expressions: 



61 - 



2[i cos(a/2) 



(-2/iVb 



. ea ea 

sm h vjt cos — 

2 2 



(12) 



2\x cos(a/2) 



ire vac 
tan — = — 
2 2S 



tt / ecu ea\ , 

Jy^y [2^Vo cos — + wt* sin —j , (13) 

(14) 



Elastic energy flow in the presence of friction a xy cx ii x 
depends |12) on the distance from the tip r: 



Ji 



fJ,C 



2„e 



(15) 



where K is a coefficient in front of the main contribution 
in the strain |9 x u x | = K\r\^ 1+£ ^ 2 near the tip in the 
slip region at y — 0. In the expression for the rear tip e is 
replaced by — e. It can be easily seen that on each scale the 
change of the elastic energy flow equals the energy loss due 
to friction. Thus, we must allow that on the microscopic 
scale a the friction law is changed so that the elastic energy 
flow on smaller scales remains constant, the same as for 
r = a. For example, this should be the case, if the frictional 
shear stress saturates at large velocities. We assume that 
this microscopic cutoff a is a material property which can 
be considered to be the same for the head and rear tips 
and essentially is independent of e. Then we can express 
the energy flows in terms of Vo, t* and e: 



•h = 



■h 



(2/iVo ' 



(l+e)a 



(l+e)ct 



(ha) 



sin(7re) 



(2 sin a) 1+e w/j,k 



(2mV c 



(l — e)a * • f 

COS — k 2 WT Sill 1 



(16) 

^~T^) ( fca ) 6 sin(Tre) 



(2 sin a) 1 e w/j,k 



We argue that if the height is larger than either of 
the lengths X s i or \ s i and smaller than the other one, 
the solutions for large height are valid essentially. Minor 
changes can be determined with the help of the solutions 
for isolated stick and slip regions given in Appendix iBl 



3.2 Large wavelengths, A >> H. 

In the limit of long slip and stick regions, one obtains 
"zero-mode" solutions with the relaxed stresses modified 
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by the presence of friction. The solutions are valid far from 
the slip edges. In the stick region |x| < X s t/2, 



V 



[H - y)x 
H 



V 



y(y - H) 

4(1 - v)H 



>xy 



(17) 
(18) 



and the normal stresses ~ /j,Vq/(1 — 2v) depend only on 
y. Here Co is related to the displacement uq of the sample 
top with respect to the middle of the stick: <to = /j,uq / H. 

In the slip region shear stress is relaxed to fiv/S and 
the sliding velocity is equal to the driving velocity v with 
exponentially small corrections far from the boundary 
points. 

Since, in the limit of small height, the normal stresses 
are relaxed in the slip region, we can find the energy fluxes 
similarly to the case of semi-infinite slip |12) . 

Using the asymptotic solutions described above, we get 
energy flows to (from) the regions of length ~ H surround- 
ing boundary points: 



3.3 Stability and degeneracy 

Solving the elastic problem we have found energy flows 
to the boundary points of slip pulses Ji i2 - Equations of 
motion for the tip and the resticking end of the crack are 
given by the energy conservation, Eqs. I)10[l. Therefore, 
we cannot determine four quantities (c, t* , k = 2tt/X, 
i] = a/n = (A — A s t)/A) using only two equations. Two of 
the parameters cannot be chosen directly, and we are faced 
with the problem of dynamical selection. We will address 
this issue in the next section using stability arguments. 

There is another kind of degeneracy related to the fact 
that the sign of the term in the parentheses in Eqs. (|16fl 
and iffi])) for J 2 can be either negative or positive, depend- 
ing on relative strength of the average shear stress and 
velocity. For simplicity, we consider this question in the 
case of large wavelengths. For small wavelengths, the con- 
clusion will be similar. In the following we neglect the tip 
dissipation (c/co << 1) assuming that the main dissipa- 
tion comes from the friction in the slip region. Then Eqs. 
(TTUj) reduce to J\ = J 2 = 7c. The velocities of the rear tip 
(ci) and the head tip (c 2 ) of the stick region are (see Eq. 
© and CT) 



(0) [icH Zero , X st 
Jl < 2 - ~2~ [j ± 2~H Vo 



(19) 



where the subscripts 1 (upper sign) and 2 (lower sign) 
correspond to the regions near the rear and head tips of 
the stick region, respectively. 

Without friction, invcrse-square-root singularities ap- 
pear in the strains and stresses within a distance ~ H 
around the boundary points. The friction changes the ex- 
ponents to — l/2±e and creates additional dissipative en- 
ergy flows from the sample. These dissipative flows are also 
restricted within the distance ~ H in the slip regions. 

Thus, the expressions for the energy flows, instead of 
Eq. Q16[l . take the following form in the limit of small 
height: 



HcH ( a X st \ ( H 
Jl = — {j + 2H V °) U 

_ ficH ( cr X st \ 2 / H 
'- — { J~2H Vo ) U 



(20) 



These equations can be rewritten in terms of the average 
shear stress t* using the relation 



_X_ 

Xst 



(21) 



Strictly speaking, if 1 — 2v <C 1, the relations derived 
here are not valid on intermediate scales H <C X st , X st 
H/ \/l — 2u. In AppendixOwe show how to treat the elas- 
tic problem on all scales larger than H without friction. 
The results on intermediate scales do not substantially dif- 
fer from those obtained above (only by numerical factors 
if written in terms of r*). 



ci 



TrSlnZii 
w ln(JJ/a) 



C2 



ttS In Ai 
vj\n{H/a) 



where 



H ( fj,VoX st 

Al > 2 = 2^{-vr ±(To 



(22) 



(23) 



For the steady-state motion c\ = c 2 = c and Vo = v/c. 
There are two branches of the solution. On the first 
branch, we have: 



2H 



— a n = <T„ 



(24) 



thus fiVoX s t/(2H) > cr c , and cto can be arbitrarily small. 
On the second branch, 



4 



2H 



(25) 



This branch contains the special static solution c = v = 0, 
Vb = for (To = <r c . This solution is unstable with respect 
to a change of A s * . The instability can be seen from the fact 
that the equilibrium corresponds to the energy maximum. 
On the contrary, for the special static solution with uq = 
0, c = on the first branch, the total energy of the stick is 
dominated for small X s t by the surface energy —jX s t, and 
for large X st by the elastic energy 



W2 



H 



dx ■ 



1 / iiVqx 



-X a t/2 



2/i 



H 



LiV 2 X 



2\3 

st 



24H 



(26) 



and hence there must be a minimum of the total energy 
corresponding to the stable equilibrium. 
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A more general stability criterion can be derived as 
follows. Consider the effect of a small change of the stick 
length SX st . We suppose that resticking occurs with the 
same strain Vb- The change of the stick length leads to the 
changes of the tip velocities: 



Sci - 
5c 2 

• : hence, 
5c2 — Sci 



2nS 



VV 



w\n(H/a) [iVoXst + 2Hao 
2irS (j,Vq 



w\n(H/a) /iVoXst — 2Hoq 



SXst , (27) 
SX st , (28) 



2tt5 



2(^o) 2 A st 



wln(H/a) IVoXst) 2 - {2Ha y 



(29) 



The necessary stability condition (Sc2 — 5ci)/5X s t < is 
satisfied only on the first branch. 



4 Dynamical selection of the stick-slip pattern 

We note that no deterministic approach can fix the two 
free parameters which appear in the result. In this section, 
we discuss the selection problem for the parameters of the 
stick-slip pattern on the example of a long sample, starting 
from the limit of small driving velocities. 

We focus our attention on the first stable branch. For 
small heights, different stick regions do not interact, and 
each of them is characterized by its length X s t , propagation 
velocity c and shear stress parameter ctq (see Eqs. (|22l24tl 'l : 



hVqX s 
2H 



v 

V = -. 

c 



2irS 
\n{H/a) 



In 



CTO 



(30) 



In the exactly stationary regime c is the same for different 
stick regions. The solutions (c, X s t) can be parametrized 
by cto. The wavelength A is a free parameter. 

It is plausible to assume that do cannot be much larger 
than a c . The inequality <jq 3> <J C would lead to large val- 
ues of the shear stress a ^S> cr c in regions of the length ex- 
ceeding H . Such a region is unstable with respect to slip. 
Larger values of Co are not impossible, they are merely 
less probable. Thus, we can find an upper estimate for the 
crack tip velocity c max : 



2nS 

Cmax ~ \n(H/a) ' 
and a relation for the stick length: 



A* 



2Hca c 

fJ,V 



< X 



st,rnax 



4nHSa c 
ym ln(H/a) 



(31) 



(32) 



We see that the upper limit for X s t vanishes with in- 
creasing v. In the case of small heights, only X s t > H is 
possible and hence v < v* , where 



4irSa r 



Anv r 



fj,]xi(H/a) bi(H/a) 



(33) 



Concerning the value of X s t, two cases can be con- 
sidered. The first possibility is to have c ~ c max and 

Xst ^ Xst,max- 

The second possibility following from Eq. (|30|) is c ~ 
wKK) < c m ax corresponding to a < &c which 
would lead to 



X s t 



Ijlv hx{H/a) 



< X s 



(34) 



If long stick regions are present in the sample, they will 
eventually catch up with shorter and slower stick regions, 
and then probably merge with them. Consequently, it is 
plausible that the lengths of the stick regions are mostly 
of order of X s t,max, specified in Eq. and r* ~ a c . As 
v decreases, this length diverges. 

On the other hand, the length X s i should be of the 
order of H for small v because large slip regions, X s i >> 
H with small frictional stress would be unstable against 
resticking. The size of the critical nucleus of the slip inside 
the stick phase with Co ~ a c is of the order of H which 
also prevents X s i to be much smaller than H. 

On the opposite side of the stick-slip regime as v ap- 
proaches v c , the same arguments lead to A S f ~ H. Thus, 
we can employ the condition 



min{A st , X s i} — H. 
This is fully compatible with the condition 



(35) 



(36) 



which we have found in the small velocity limit and also 
using the analogy to the two-phase region in the first order 
phase transitions. 

The results shown in Figures - represent numer- 
ical solutions of the equations of motion, J\ = 7c and 
Ji = 7c, together with two additional constraints I|35I36() . 
We model the first condition by the requirement that 
kH/siaa is a fixed number. The dependences of the slip 
fraction and the dimensionless crack tip velocities on the 
driving velocity are evaluated both for the limits of large 
and small heights. In Fig. |S] the values of \ st and A are 
given in units of H . The crack velocities c remain of the 
order of S for all values ofv/v Cl the fraction of slip r\ grad- 
ually changes from to 1 and the wavelength A diverges 
at the both ends of the stick-slip region being still of the 
order of H for the intermediate values of v/v c . 

The sliding velocity in the periodic stick-slip regime 
calculated from Eqs. (|11I13|) is presented in Fig. ©for a 
certain value of the driving velocity. Note that the singu- 
larities at the head and the rear of the slip pulses are quite 
different. 



5 Discussion 

A low velocity dynamics consisting of periodic slip pulses 
which heal at a critical value of the local slip velocity is 
usually attributed to the existence of a V-weakening fric- 
tional regime 5 . In this unstable regime the sliding ve- 
locity decreases with the increase of the stress. Our model 
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0.1 - 



1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

v/v e 

Fig. 3. Ratio of the crack propagation velocity to the velocity 
scale S set by viscous friction as a function of v/v c for two 
values of kH/(jr sin a): 2 (solid line) and 0.4 (dashed line); For 
all plots ira/(2H) was 2 • 10~ 6 . 



1 




0.2 0.4 0.6 0.8 1 

v/v c 

Fig. 4. Fraction of the slip ij = h/t as a function of v/v c 
for two values of kH/(n sin a): 2 (solid line) and 0.4 (dashed 
line); rj does not strongly depend on kH / (ir sin a) . For all plots 
Tva/(2H) was 2 ■ 10~ 6 . 

does not contain intrinsic instabilities of the friction law as 
a V-weakening frictional regime. The homogeneous slip is 
still metastable below the critical velocity in our picture. 
Critical shear stress and critical velocity depend not only 
on material properties but also on the geometry and they 
are related to the point of discontinuous transition (to the 
Griffith threshold for crack propagation). These two sce- 
narios do not contradict each other and can be combined 
in one model (see, for example, |14| for a description of 
"shear induced melting"). This is the usual situation in 
first order phase transitions: both, the transition point 
and instability points exist, forming the metastable re- 
gions around the transition point. The characteristic ve- 
locity scale for the discontinuous transition and velocity 
scale for a V-weakening unstable frictional regime can be 
well separated leading to stick-slip even above the insta- 
bility point but, of course, below the transition velocity. In 




1 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 

v/v c 

Fig. 5. \/H (solid line) and \ s t/H (dashed line) as functions of 
v/v c for kH/(nsina) = 0.4. For all plots na/(2H) was 2-10" 6 . 



■s 




0.4 - 
0.2 - 

I ' 1 1 

-1 -0.5 0.5 1 

x/X 

Fig. 6. Dimensionless sliding velocity U;c|j, = o in the peri- 
odic stick-slip regime as a function of the coordinate for 
kH / (tt sin a) = 2. Two periods are shown. The sliding velocity 
is zero in the stick regions around |jc| = A/2. Pattern parame- 
ters are calculated for the driving velocity v = 0.5v c . 

our model we discuss the limiting situation where the V- 
strcngthcning regime holds even for very small velocities 
(at least smaller than our transition velocity. 

Now let us discuss the relation of the obtained re- 
sults to the very interesting experimental observations of 
Baumberger, Caroli and Ronsin They performed ex- 
periments of a gel sliding on a glass plate. The driving 
velocity was given and the shear force, and thus the av- 
erage shear stress, was deduced from the spring elonga- 
tion. Above some critical driving velocity v c ~ 100/im/s, 
steady sliding was observed. At velocities smaller than the 
critical one, periodic stick-slip sets in (see figures in 
IS])- Upon decreasing the driving velocity v, no hystere- 
sis of the transition was detected. In the stick-slip regime 
they observed the propagation of self-healing pulses with 
no opening, nucleated periodically at the trailing edge of 
the sample. The propagation velocity of these cracks was 
about 60 times larger than the critical sliding velocity, yet 
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still much smaller than the shear wave speed. The crack- 
like singularities in the slip velocity field were detected 
behind the tip of the pulse but not at the resticking end. 

The authors very carefully studied the dependence of 
the critical velocity, crack velocities and other properties 
of slip pulses on the properties of the gel, but the depen- 
dence on geometry was not discussed. Since both the criti- 
cal velocity and the critical stress, predicted by our model, 
depend on the height of the sample and not only on ma- 
terial parameters, the crucial test for the relevance of our 
model to this experiment is not possible yet. Moreover, we 
assumed a very large sample length compared to its height 
and discuss a spatially periodic stick-slip pattern, while in 
the experiment at most one slip pulse was simultaneously 
observed in the sample nucleated at the trailing edge. The 
length of the slip pulse increased with the driving veloc- 
ity v and sometimes the sample was too short to observe 
more than one stick-slip boundary. Furthermore, the fact, 
that no hysteresis of the transition was detected in the ex- 
periment, points towards the possible interference of the 
intrinsic instabilities for low sliding velocities, which have 
not been discussed in our model so far. However, as we 
have already mentioned, the characteristic velocity scale 
for the discontinuous transition and the velocity scale for 
the V-weakening unstable frictional regime can be made 
well separated, if smaller heights of the sample would be 
available for experimental investigation. 

The observed nonlinear behavior of the stress with 
the velocity for relatively high sliding rates (the so-called 
shear-thinning rheology) can, in principle, be incorporated 
into a more sophisticated version of the theory presented 
here. This has already been discussed shortly in the con- 
text of the microscopic cutoff a, the small region around 
the crack edge where the transition from stick to slip takes 
place. 

Another way to approach the observed dependence of 
the shear stress on the sliding velocity is to make a simple 
modification of the friction law in the slip region which 
amounts to the addition of the constant: a xy = const + 
/3ii x . This modification does not lead to any changes in 
our results except for a trivial shift of all values of the 
shear stresses. 

The predictions of our theory may differ from experi- 
mental results due to different geometries, interference of 
the intrinsic instabilities, some specific features caused by 
polymeric nature of the gel etc. The effects on the prop- 
agation of single slip pulses are planned to be considered 
elsewhere [*T7] . 

Our aim was to construct a generic, conceptually sim- 
ple model which allows for stick-slip motion in sliding fric- 
tion without call to some very specific properties of mate- 
rials. In our model all complicated properties of the gel are 
hidden in only two material parameters which are related 
to surface properties: the energy scale 7 and the velocity 
scale S due to the viscous frictional law. Of course, these 
parameters may depend on the applied normal stress but 
we assume that they exist also for zero normal stress. The 
velocity scale S can be easily estimated from the shear 
stress response to a jump of the driving velocity in the 



sliding regime (see Fig. 15 in 5 ). This response is pure 
linear in our model and leads to 

a(t) = a fi n - {a fin - a in ) exp(-St/H), (37) 

where <7j„ and cr/m are the initial and final stresses be- 
fore and after the jump. The response is not changed if 
one takes a modified friction law by adding a constant as 
described above. This exponential behavior is in a good 
agreement with the experimental curve for large changes 
of stresses and velocities (see Fig. EJ) and gives an esti- 
mate for S/H w 1 s _1 , which leads to S w 1 cm/s for 
used samples with H — 1 cm. 



0.8 




0.1 ' 1 1 1 1 1 1 

1 2 3 4 5 6 

t(s) 

Fig. 7. Several experimental points taken from Fig. 15 of Ref. 
5 : relaxation of shear stress after a jump of the driving veloc- 
ity. The line represents a fit according to Eq. 1371 . 

One would expect that for ordinary elastic materials 
the velocity S should be of the order of the speed of sound. 
However, for gels the shear modulus /1 is much smaller 
than for ordinary materials. The shear sound velocity c t = 
(fj,/ p) 1 / 2 is only 2 m/s. The velocity S = fi//3 is linear in [i 
and should be even smaller. This is a possible explanation 
for a relatively small value of S compared to Ct- 

The existence of a critical sliding velocity, where the 
transition to a stationary sliding occurs, appears natu- 
rally in our description and is given by Eqs. II) . The 
length of the slip pulses increase as v approaches v c , and 
the fraction of the slip phase tends to unity in qualitative 
agreement with the experiments. The characteristic value 
of the shear strain in the sliding mode near the critical 
velocity experimentally was about u xy — 0.04. Thus, we 
can estimate from Eq. the characteristic difference be- 
tween the interface energies in the slip and stick states to 
be 7 = 0.1 J/m 2 . 

In the stick-slip regime, which exists below v c , the nu- 
cleation of a slip pulse takes place at the trailing edge of 
the sample and requires overshooting above the Griffith 
threshold according to the experiment. This overshoot- 
ing is not small and we can use Eq. JSJ to estimate the 
crack velocity. Because of the week logarithmic parameter 
dependence, we conclude that the crack propagation ve- 
locity c is of the order of S and essentially independent of 
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the driving velocity, in agreement with experimental ob- 
servations. After the slip pulse traversed the sample the 
stress drops below the Griffith threshold and resticking 
takes place via propagation of a healing front. Theoretical 
predictions for the crack velocities c in the steady-state 
spatially periodic stick-slip regime are given in Fig. |3J we 
get c ~ S ~ 1 cm/s, which is also in qualitative agreement 
with the experiments. 

The crack-like singularities in the slip velocity field are 
predicted by our theory and were detected behind the tip 
of the pulse. At the resticking end the sliding velocity in 
the experiment suddenly drops to zero from a finite value 
which is close to the critical velocity. We should note that 
this discontinuity with finite amplitude is not compati- 
ble with solutions of the elastic problem which allow ei- 
ther crack-like singularities or, as a particular solution, 
a smooth, continuous termination of the slipping process 
at the healing point. This type of solutions, which corre- 
spond to a vanishing stress intensity factor due to some 
additional constraint on the parameters, was discussed in 
the literature (e.g. Refs. ^HIE]) and also in the previous 
sections. Perhaps, one could argue that a fast transient, 
too fast to be experimentally resolved, takes place (see 
Fig. El- 

The geometry of the experiment is such that the total 
macroscopic friction of the sliding sample and the nucle- 
ation of the pulses depend on the processes taking place 
at the edges of the sample. The stresses here are highly 
inhomogeneous and the kinetic phenomena should be con- 
sidered with a great care. It would be very interesting to 
perform a experiment in such conditions where the length 
of the sample and its height are well separated and to 
study the possible height dependence of the critical veloc- 
ity and the shear stress. We believe that further theoreti- 
cal and experimental investigations will shed light on this 
phenomenon where two intriguing problems, crack propa- 
gation and friction, combine together. 
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We consider only the plane strain case with u z = 0. 
The solutions can be found in form of Fourier series jll) : 



=vt + u*{y) +ReJ2 (a m e- ms - v + b m e- ms + y 

771— 1 

+c m e ms - y + d m e ms + y ) e™^"^), 



(39) 



m—l 

ik 



ik 

&m 6 



-b e 



c m e ms - y - l -^±d m e ms + y ) e imk ^- ct \ (40) 
k 



&xy =t* + A^Re mk 



m—l 



s_ k 



-1-2— d e ms+y ) e imk ( x - ct ) 
k m J 



(41) 



where 




and k is the pattern wave number: k = 2n/ X. In a periodic 
pattern, an average value of a quantity corresponds to its 
zero Fourier harmonic. We assume that on average there is 
no normal stress, i.e. u*(y) is a linear function, and there 
is no zero-harmonic term in u y . 

We use boundary conditions u y \ y =o — u v \ v —h = 
and d x u x \ y =H = 0, which give three linear uniform equa- 
tions for a m , b m , c m , d m . This allows to express the stress 
t = <J X y\ y= Q in terms of the Fourier harmonics of the di- 
mensionless sliding velocity V = u x /c = Vq — d x u x : 



T{rj) =t*+ ^Rc J2 (-2i)x(m) J B r) 



G imk( X — ct) 



m—l 



where B m is 



1 



TT/k 



B m = — I dxkV(x)e 

7T J 

-n/k 



A Relation between strain and stress at y = 

We consider only periodic stationary solutions, possible 
generalizations are specially mentioned. 

A starting point is the momentum conservation equa- 
tion of the linear elasticity theory: 



d t u a 



(cf -c 2 t )V a Vf3U(3 + c 2 t V 2 u a , (38) 

where c/ and Ct are the longitudinal and transverse sound 
velocities respectively, c;/c t = ^/2(1 — v)/(l — 2v). 



and x( m ) is calculated under the assumption c 2 <g; c|, cf: 



X = 2(l-i/) 



(3 — 4z/) sinh(mfc_ff ) cosh(mfciJ) — mkH 
(3 - Av) 2 sinh 2 {mkH) - {mkH) 2 



This is the result of the static elasticity theory. The first 
c-dependent correction would be proportional to c 2 /c 2 . 

We note that the limit of v — > 1/2 and mkH — > is 
not well defined: the limiting value of \ depends on the 
ratio {mkH) 2 /(l - 2v). If {mkH) 2 < (1 - 2v) one gets 
x(m) = l/{2mkH). The relation between zero harmonics 
of the shear strain and the stress r = 2/iu xy is recovered 
on scales larger than Hj\J\ — 2v. If v is close to 1/2, there 
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is an intermediate region y/1 — 2v -C mkH -C 1 where we 
have x(m) « 2 /(mkH). Finally, if mkH ^> 1, 



2(1-^) 
3-4^ 



The latter is realized for the periodic pattern with the 
wavelength smaller than the height, and all the terms in 
the Fourier series satisfy the condition mkH 3> 1. 



B.2 Single stick of length X st < H 

The solution with finite friction reads: 



bix-\-b2 



V(x) = { (x+A st /2) — (x-X st /2) 



jjr, |x| > X st / 2 

2 

0, \x\ < X st /2 



Here 



B Short isolated stick and slip regions 

In this Appendix we specify the results for short isolated 
slip and stick regions. This limit is characterized by the 
following relation between the strain (expressed in terms 
of dimensionless velocity V in the stationary regime) and 
the stress t(x) = cr X y\y=o'- 



7TH7 



dz 



V(z)+r\ 



(42) 



which is derived similarly to the periodic case. 



B.l Single slip of length A s ; <C H 

First consider the case of a short slip, propagating with 
speed c, in the presence of friction proportional to sliding 
velocity; the top surface is uniformly displaced and fixed 
(no driving velocity). 

The solution with zero friction which satisfies the con- 
dition V = in the stick region |x| > X s i/2 has the fol- 
lowing form: 



bix+b2 



M < Xm/2 



V(x) = < V(W2) 2 -* 2 ' (43) 
1 0, \x\ > X sl /2 



where 



VJT 



1 



(u 2 - Ml) , 



2\l ' 7T 
U\ = «;eU=A s ,/2 , U 2 = U x \ x= _x sl /2 ■ 

The solution with the linear friction law reads: 



bi x+b2 



V(x) 
where 



(\ sl /2-x)— (\ Bl /2+x) — 



T±r, |x| < X s i/2 
0, \x\ > Xsi/2 



tan 
b 2 = 



7T£ WC 

~2 ~ 2S' 

ewT*X s i cos \ 



2(i. 



„ COS -5- 

— + {U2~Ul) =- 



4/i ' 7T 

Extracting the coefficients in front of singularities, one 
finds the energy flows: 

2 



c (2fi(u 2 - «i) + ™ T 'Y e±1) ) S in(^ e ) ( a x 1 ' 



A-K 2 VJ[lXs 



X s 



The energy flows could have been recovered from ex- 
pressions Ultty for periodic stick-slip in the limit a — * 0, 
2a/k = X s i = const. 



7re vjc , wr* cot ^ 

tan — = — , b\ = Vq = — 

2 2S*' 2^ 



b, = — 



weX.tT* cot %r 



4/i 



(^x I x — + oo — oo) 



cot^ 



C Small height solutions for almost 
incompressible media, 1 — 2v 1 

The basic idea of the small height approximation is to use 
the approximate equations of elasticity theory which are 
valid if the transverse gradients (along the shortest verti- 
cal dimension) of displacements are much larger than the 
longitudinal ones. This is justified by the fact that certain 
boundary conditions are imposed on the upper and lower 
surfaces. The approximation is adequate in the regions 
far from the points of discontinuity in the boundary con- 
ditions. In those regions, imposing u y on the boundaries 
additionally implies that d x u y -C d v u x . The approxima- 
tion cannot be used at the distances of the order of H near 
the crack tips. 

Neglecting d 2 terms in comparison with dy terms in 
the exact equations, one gets: 



\2d x (d x u x + dyUy) + bldyU x = 

!2dy(d X U x + OyUy) + l^^Uy = , 



(44) 
(45) 



where &q = 12(1 — 2v). The case of bo ~ 1 was discussed 
in Sec. 13.21 In the limit bo <C 1 one can conclude from Eq. 
(|45fl that divit = d x u x + d y u y is independent of y. Then, 
it is easy to see that u x quadratically depends on y, and 
u y is a cubic polynomial in y. 



C.l Long stick 

In the stick region, we use stationary boundary conditions 
to find the solution which is expressed in terms of unknown 
functions of x: 



divu = fx(x) 

u x = y(y - l)gi(x) + u y - V x(y 
u v = v(v - i)(pi(aO + yqi(x)) ■ 



i) 



(46) 
(47) 
(48) 



Here the lengths are measured in the units of height, and 
Uo is the displacement of the upper boundary with re- 
spect to the origin. Together with the first equation of 
(|45fl . there are four linear differential equations for the 
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four functions, which result in g'{ — b^gi — 0. Thus, all 
functions can be parametrized by two constants: 

gi = asmh(b x) + bcosh(b Q x), /i = -Vq — -g[ , (49) 
qi = -~g[, Pi = V - fi = ~V Q + ^g[ . 



C.2 Long slip 

The same procedure can be used for the solution in the slip 
region. Here, because of different boundary conditions, u x 
has another form: 

u x = u + ^V (y -l) + (y- l)(h 2 (x) + yg 2 {x)) . (50) 

The presence of the boundary friction (finite S) makes this 
problem more difficult. The set of equations is reduced to 
the third-order differential equation: 



acosh(6 A s t/2) - 2Acosh(& A ;s z/4) 



4h';-b 2 h 2 + -(-h 2 " + b 2 h' 2 )=0. 



(51) 



Without friction the solution can be written as 

h 2 = Asmh(b x/2) + B cosh(b x/2) . (52) 

Small friction changes the balance between the growing 
and vanishing exponentials on large scales but the final 
effect can be shown to be negligible. 

Functions g 2 , p 2 , <li and f 2 = divu are expressed in 
terms of h 2 , Vq and c/S, in analogy to Eqs. (@H1>- 

C.3 Matching stick and slip solutions 

The solutions found above are good approximations far 
from the stick-slip boundaries. To match the solutions, one 
must satisfy global mass and momentum conservations. At 
every boundary point this will give two relations, which 
equals the number of unknown constants in every interval. 

The condition, that the total horizontal force acting on 
a piece including the stick-slip boundary is zero, leads to 
the continuity of divergences: /i(±A st /2) = f 2 (T^si/2)- 

The second condition is the equality of the volume in- 
tegral of divergence to the surface integral of the displace- 
ment. This gives the relations 

±W X st - i<?i(±A st /2) + ~h 2 {T\ s l/2) = \u . 
4 o 6 2 

Upper and lower signs correspond to the left and right tips 
of the slip region. 

Deriving the conditions, one neglects 1/X s t and bo in 
comparison with unity. The full set of equations for four 
constants reads: 



asinh(6 A st /2) + AA sinh(& A s z/4) 



3V \ S 



bo 



6sinh(6 A st /2) + 2B smh(6 A si /4) =0. 
Energy flow to the rear edge of the slip is given by 



M c ( 4 



1, 



J= ^l l B z -^b z -ut- 



2Voa{cos,h(boX st /2) - 1) 
bo 



2V bsmh(b X st /2) 1 2 2 
7 + Vou X st - -Vq X st 

D 4 

In the limit of not extremely long stick and slip regions, 
H <C X s t, X s i <C H/bo, the expressions for the energy flow 
are reduced to 

^ C fir ^ , 2 ( A s« + 

' 1 1 ■ ■ . (53) 



^ = M[ VoX * t±U °-4X st+ X sl 



Relating the displacement uq and the average shear stress 
t* as 

4A st + A s( Ht* 
4A st fi 

we rewrite the energy flows in the following manner: 
ficH (VoX st t* X 



6cosh(6 A s t/2) — 4_Bcosh(fe A s ;/4) = — 3u , 



M 2A st/ 

This result differs from the frictionlcss case of Eqs. (|20|l 
only by numerical factors. Moreover, for A = X s t and Vo — 
0, the energy conservation, J = 7c, leads to the result for 
the critical stress: r* = 2a r . 
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